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0> ■ Abstract 
O , 

^ ■ In the framework of the gluon-gluon fusion process for Higgs boson pro- 

Xf^ . duction there are two different prescriptions. They are given by the exact 

^ ! process where the gluons couple via top-quark loops to the Higgs boson and 

i-C \ by the approximation where the top-quark mass rut is taken to infinity. In 

the latter case the coupling of the gluons to the Higgs boson is described by 
an effective Lagrangian. Both prescriptions have been used for the 2 — > 2 
^ . body reactions to make predictions for Higgs boson production at hadron col- 

liders. In next-to-leading order only the effective Lagrangian approach has 
^ ; been used to compute the single particle inclusive distributions. The exact 

;h ' computation of the latter has not been done yet because the n-dimensional 

^ extensions of 2 — > 3 processes are not calculated and the two-loop virtual 

corrections are still missing. To remedy this we replace wherever possible 
the Born cross sections in the asymptotic top-quark mass limit by their ex- 
act analogues. These cross sections appear in the soft and virtual gluon 
contributions to the next-to-leading order distributions. This approximation 
is inspired by the fact that soft-plus-virtual gluons constitute the bulk of 
the higher order correction. Deviations from the asymptotic top-quark mass 
limit are discussed. 
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1 Introduction 



In the past few years many articles have appeared on searches for the Higgs 
boson and the reactions in which they are produced. One of them is the 
gluon-gluon fusion process. According to the standard model gluons do not 
interact directly with the Higgs boson but the coupling is mediated by a 
fermion loop. Since the coupling of the Higgs boson to fermions is propor- 
tional to the mass of the fermion the reaction proceeds mainly via a top-quark 
loop [1]. The lowest order loop is a triangle graph and the Higgs boson decay 
rates into two gluons or two photons were already calculated at the end of 
the seventies [2]. The first calculation in the gluon-gluon- fusion model for 
the production process was done at the end of the eighties by [3]- [5], (see 
[6], [7] for later references). Reactions hke g + g^g + H,q + q^g + H and 
q + g ^ q + H were calculated. In particular the first reaction involves a box 
diagram leading to complicated dilogarithms already on the Born level. In 
the early nineties people succeeded in calculating the next-to-leading order 
(NLO) corrections to the total cross section which involved the computa- 
tion of the two- loop triangular graph with an external Higgs boson [8] . The 
calculation could be greatly simplified by taking the infinite top-quark mass 
limit. In this limit the gluons couple directly to the Higgs boson and the 
Feynman rules are given by an effective Lagrangian. It turned out that the 
latter method gives a good description of the exact calculation [9] provided 
the Higgs boson mass niH and the transverse momentum pt are smaller than 
the top-quark mass rrit [4], [5], [7]. In particularly the total cross section 
receives its main contribution from small pt- If the Higgs mass is not too 
large {mn < 2 rrit) the effective Lagrangian gives a good description of the 
total cross section so that recently one has also finished the next-to-next- 
leading order (NNLO) computation [10]-[14]. However at Higgs masses and 
transverse momenta equal or larger than the top-quark mass the differential 
cross sections calculated with the effective Lagrangian method start to devi- 
ate from the exact cross sections. This has been checked on the Born level 
in [4] , [5] , [7] . The investigation should now be done in NLO but we realize 
that the exact cross sections are not available. Differential distributions in 
NLO using the effective Lagrangian (or the — > oo approach) have been 
calculated in [15]-[18]. In the same approach the resummation of the loga- 
rithmically enhanced contributions to da/dpt at small pt have been carried 
out in [19] -[21]. The first landmark calculation to get the full NLO differ- 
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cntial distribution has been achieved in [22]. In the latter one has exactly 
calculated all matrix elements of the 2 — 3 processes. These reactions even 
contain one-loop five-point functions. However the calculation of the graphs 
uses the helicity method in four dimensions. To compute the single particle 
inclusive process we need the matrix elements in n dimensions. Moreover 
the two-loop virtual corrections, which arc needed to cancel the infrared and 
coUinear divergences, have not been calculated yet. Therefore we propose 
to make an approximation by replacing all Born contributions in the infinite 
top- quark mass limit by their exact analogues in the virtual-plus-soft correc- 
tions. However this is not sufficient. We have also to demonstrate that the 
soft-plus- virtual gluon approximation gives a good description of the differ- 
ential cross section. Using a certain prescription we can show that this is 
really the case. 

Our paper is organized as follows. In Section 2 we present the formulas 
for the exact cross sections and their analogues in the infinite top-quark 
mass limit. Then we make approximations for the partonic soft-plus- virtual 
and the soft-gluon cross sections. Finally we adopt a prescription how to 
implement these formulae for the hadronic pt distributions. In Section 3 
we make comparisons between our approximate differential distributions and 
those which are derived in the limit mt — > oo. 
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2 Approximation to the exact differential cross 
section for Higgs production 

The differential process we study is the semi-inclusive reaction with one Higgs 
boson H in the final state 

ifi(Pi) + H2{P2) ^ H{q) +' X' . (2.1) 

Here Hi and H2 denote the incoming hadrons and X represents an inclusive 
hadronic state. In our study we limit ourselves to 2 ^ 2 and 2^3 partonic 
subprocesses. The kinematics of the 2 — > 2 reaction is 




Figure 1: The exact process g + g ^ g + H. 



a{pi) + b{p2) ^ c(p3) + H{q) , 

s = {pi+ P2f , t={pi-qf , u= {p2 - qf 



(2.2) 
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The exact calculations of the 2 — > 2 processes are given in [3]- [7]. They 
consist of the following parton subprocesses 

g + 9^g + H, q + q^g + H, q{q) + g ^ q{q) + H . (2.3) 

g + H subprocess in Fig. 1 is equal 



The Born cross section for the g + g 
to 



j2 (l),exact 

dt du 



H 



N 



16n stu 



\A2{s,t,u)\^ + \A2iu,s,t)\' 



with 



+ \A2{t,U, s)\'^ + \A4{s,t,u) 



TT 



S{s + t + u-mjj) , (2.4) 



(2.5) 



in sin^ 9w 

where e denotes the electric charge and 9w is the weak angle. The constants 
Mw and G p denote the mass of the W and the Fermi constant respectively. 
Further we want to mention that N — 3 for QCD. The dimensionless func- 
tions A2{s,t,u) and A4(s,t,u) are given in the Appendix of [4]. The Born 
cross section for the q + q ^ g + H subprocess in Fig. 2 equals 



j2 (l),exact 

dt du 



3 



12877 s{u + t)'^ m 



A5{s,t,u) 



xd{s + t + u- ttih) , 



(2.6) 
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Figure 3: The approximate process g + g ^ g + H . 



where the function ^5(5,^,^) is given in the Appendix of [4]. Finally the 
Born cross section for the q{q) + g ^ q{q) + H reaction becomes (see Fig. 2) 

72 (1), exact 9 9 4-1 

^ dtdu ~ 1287r + M2. 7vl^'^''^'''^l 

x5{s + t + u-m]j) . (2.7) 

In the limit of infinite top mass rrit the functions above simplify enormously. 
Actually they can be derived from the effective Lagrangian 

A// = G$(x)0(x), with 0{x)^-\gI^{x)G^'^^{x), (2.8) 

where $(x) represents the Higgs field and G is an effective coupling constant 
given by 

The quantity C is the coefficient function which describes all QCD corrections 
to the top-quark loops in the limit — > cxo. For external gluons, which are 
on-shell, the latter quantity has been computed up to order as in [8], [9] 
and up to in [24], [25]. Up to second order it reads 



2^ 



2 r 



2777 ii: 
h 19 In — 

18 m- 



2 
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Figure 4: The approximate processes q + q ^ g + H and q{q) + g ^ q{q) + H. 



/ 67 16 ^2 ^ 



(2.10) 



Here Hr represents the renormahzation scale and Uf denotes the number of 
hght flavours. Moreover a^^^ is presented in a flve-flavour-number scheme. 
In the infinite top-quark mass hmit the Feynman rules can be derived from 
Eq. (2.8). In that limit the Born cross sections become 



(l),mt^oo 
gg^g H 

dt du 



144 TT 



N 



1 stuM^ 



4 I J.4 I 4 I 

s +t +u +m 



H 



xS{s + t + u- rriff) , 



(2.11) 



(l),mt— >oo 
qq^g H 

dt du 



a,„ a? _ 1 ^2 



u 



288 TT 



5{s + t + u-mH) , (2.12) 



d^ 



a 



.(i).m, 
<iy- 
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dt du 



Uw ttf 1 U 



2887r N tM^ 



■6{s + t + u 



m 



H) 



(2.13) 



where the the graphs are shown in Fig. 3 and 4. The next order gluonic 
corrections to the 2^2 reactions in Figs. 1 and 2 have not been calculated 
yet. Some of the graphs are shown in Fig. 5, which shows that the calculation 
will be very tedious. However we can make an approximation. In the infinite 
top-quark mass limit the soft-plus-virtual {S + V) cross sections could be 
written as [14] 





Figure 5: Samples of two-loop graphs contributing to g + g^g + H. 



j2 " ^ab->-cH 

dt du 



^^l\[S,t,U,lX,fl )S ^^^^ 

+7r6{s + t + u-ml)^K \MB, 



/I) I 



(2.14) 



Here d'^a^^^ denote the Born cross sections in Eqs. (2.11) -(2.13) and MB^^_^^jj 
is a left over piece which is numerically very small. The term N{s, t, u, A, /j,^) 
is an universal function which depends on the parameter A which serves 
as a momentum cut off for the infrared divergence. Finally K denotes a 
combination of colour factors which vanishes in the supersymmetric limit 
Ca — Cf — rif — N. Here Ca, Cp are the standard colour factors in SU (N). 
For more details see Eqs. (5.24)-(5.26) in [14]. Since N{s,t,u, A, n"^) is uni- 
versal we replace the Born cross sections in the first term of Eq. (2.14) by the 
exact ones in Eqs. (2.4), (2.6), (2.7). In this way we get a better soft-plus- 
virtual gluon approximation for the Higgs boson cross section which is also 
valid for Higgs masses and transverse momenta pt larger than the top-quark 
mass nit. The 2 — > 3 reactions are denoted by 

a(pi) + b{p2) c{ps) + d{p4) + H{q) , 

s = (pl+P2)^ ^=(pl-g)^ ^i=(p2-g)^ S4 = (P3+P4)^ 



S4 = S + t + U — mrr . 



(2.15) 



The matrix elements for the 2 — 3 processes have been exactly calculated 
in [22] although in four dimensions. Some of the graphs are shown in Fig. 
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Figure 6: Samples of graphs contributing to g + g^g + g + H. 



6. However we need them in n dimensions to regularize the infrared and 
coUinear divergences (for n — A there is a problem see [26]). Furthermore 
wc also need the exact virtual corrections to cancel the infrared divergences. 
Since the latter arc not calculated yet we can only make an approximation 
for the soft parts (54 — > 0) of the 2 — > 3 processes. In the — > 00 limit 
these parts are 

,2 (2),SOFT , ,2 (1) 

where cPa^^'' are the Born cross sections in the limit rrit ^ 00 given in Eqs. 
(2.11)-(2.13). The term I{s.t.u.S4,fi'^) is a universal factor and contains 
simple functions which are proportional to lns4/yU^. For more details see 
Eqs. (5.16)-(5.20) in [14]. We get a better approximation to the exact cross 
sections if we replace in Eq. (2.16) the cross sections in ^ 00 limit by the 
exact ones in Eqs. (2.4), (2. 6), (2. 7). The most optimal next-to-leading order 
(NLO) cross section that one can achieve is to use the exact lowest order 
cross sections in Eq. (2.4), (2.6), (2.7) and in next order to substitute them 
in Eq. (2.14) and (2.16). This relies upon the fact that the soft-plus- virtual 
gluon approximation is a very good substitute for the exact cross section. 
We know from our experience with the cross section in the infinite top-quark 
mass limit that this is really the case. This is revealed by a study of the 
transverse momentum pt and the rapidity y distributions in Figs. 13-15 of 
[14]. Above Pt = 100 GeV/c and mn > 100 GeV/c^ the soft-plus- virtual 
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gluon approximation accounts for 80 % of the cross section. However we can 
do even better. This becomes clear if we look at the transverse momentum 
distribution 



a Pt a,b=q,g -^^-i" " 



(2.17) 



with 



_ mjj + 2p^,+2^pUpl + ml) _ 

and denotes the momentum fraction luminosity defined by 

<l«^^^(a;,/.2) = £dx,J^'dx2 6{x-x,x,)f^^{xu^i')f^\x2,^i'). (2.19) 
However Eq. (2.17) can also be cast in the form (see [27]) 

— {S,plml) = Yl ^min / dx^^r {x,iJ,^) 

" Pt a,b=q,g ^^min a;inin 



X— [X S,Pt,mH:fl ) , 

a Pt 



(2.20) 



where is the parton luminosity given by 

^T'{x,i^')^x-'^T'{x,i,'). (2.21) 

If we consider the whole cross section it makes no difference which definition 
we arc using. However if we limit ourselves to the soft-plus-virtual gluon 
approximation and moreover we set x/x^-^^i^ = 1 in Eq. (2.20) we get a differ- 
ence. In fact wc enhance the small x region which leads to an improvement 
of the approximation. This is mainly due to the fact that the small x gluons 
dominate the differential distributions as they already did in the total cross 
section (see [11]- [14]). This will be shown in the next section. 
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3 Differential distributions for the LHC and 
the TEVATRON 

In this section the hadronic differential distributions are presented for ar- 
bitrary Higgs mass ttih and top-quark mass m^. We compare the results 
for the NLO differential cross sections in the infinite top-quark mass limit 
and in the approximation derived in the previous section, which is valid for 
arbitrary m^. Since the latter is only defined for the transverse momentum 
we will limit ourselves to the pt-distributions. In this paper we will study 
Higgs boson production in proton-proton coUisions at LHC {y/S — 14.0 TeV) 
and proton-anti-proton collisions at the TEVATRON {^/S = 2.0 TeV). The 
hadronic cross section is obtained from the partonic cross section as follows 



S 



(I 1 a U a,b=q,g "'^l.min ''"1 "'^2,min ^2 

(P aab 



xff^{x2, /^^) ^ ^ ^J ^s, t, u, mfj, f/) . 

(3.1) 

In analogy to Eq. (2.2) the hadronic kinematical variables are defined by 

S^(P, + P,)\ T^{P,-q)\ U^{P2-q)\ (3.2) 

where Pi and P2 denote the momenta of hadrons Hi and H2 respectively (see 
Eq. (2.1)). In the case parton pi emerges from hadron Hi{Pi) and parton 
P2 emerges from hadron 7^2(^2) we can establish the following relations 

Pi = Xi Pi , P2 = X2 P2 , 

s — xi X2 S , t = xiiT — rn^jj) + m\ , u = X2{U — mjj) + , 

—U —Xi{T — m'^jj)—m% 

^1.-- = S + T-mj, ' " XiS+U-ml ' ^ ' 

Prom Eq. (3.1) one can obtain the pt and y distributions. Neglecting the 
masses of the incoming hadrons we have the following relations 



T — rn^j — \/S yPt + rn^j cosh y + \/~S ypj + rn\j sinh y , 
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U = rn^jj — \fS ^JpI + m?jj cosh y — ^/S ^^pf + mjj sinh y , (3.4) 
so that the cross section becomes 

The kinematical boundaries are 

C 2 

m^-^<r<0, -S -T + m]j <U < ^ + m]j , (3.6) 

T-mjj 

from which one can derive 



"^'^ ^^'-= 45cosh^y ""^"' ^'-'^ 



or 



2 (-S - mjjf 2 

-l/ma^ < 1/ < 2/max , < < ^-^ = Pj.^^^ , 

1 l + y/T^ 4 ^ (pg + m%) 

With W=2ln3--^_=, (^ + ^2 • (3-8) 

We can perform the integral over the rapidity and obtain the transverse 
momentum distribution 



yS,Pt,mH)= dy- —{S,p^,y,mH), (3.9) 

dpt dy 



with i/max given in Eq. (3.8). An alternative way to obtain the distribution 
above is given in Eq. (2.17). We checked that both procedures lead to the 
same numerical result. 

We define what we mean by leading order (LO) and next-to-leading order 
(NLO). In the infinite top-quark mass limit and in the exact computation 
the differential cross section in LO is defined by 

"-^{S^plml) = "^{S^plml) , (3.10) 
dpt dpt 
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where we shall denote the LO cross section in the infinite top-quark mass 
limit by da^'^'™*'^°° / dpt- The partonic cross sections in the latter quantity 
are given in Eqs. (2.11)-(2.13). The exact LO cross section is represented 
by da^'^'^'^^^ / dpt with the partonic cross sections in Eqs. (2.4), (2.6) and 
(2.7). The gluon-gluon-Higgs coupling is given by G in Eq. (2.9) with C = 1. 
The top-quark mass is given by nit = 174.3 GcV/c^ and the Fermi constant 
Gf = 1.16639 X lO-^GeV-^ = 4541.68 pb in Eq. (2.5). We also adopt the 
leading logarithmic representation for the running coupling and the parton 
densities. For the latter we choose the parametrizatiou according to [28] 
(namely set lo2002.dat) with A^^ = 167 MeV and nj = 5. 

The NLO corrected difi^erential cross section in the asymptotic top-quark 
mass limit is given by 



d a 



NLO.mt^oo 



dpt 



-{S,Pt,m]j) 



1 + 22 ' ^ ' 



d, a 



LOjiit^oo 



dpt 



-{S,p1,m%) 



+ 



dpt 



-{S,pl,m]j) 



(3.11) 



In d all partonic cross sections use the asymptotic top-quark mass 

limit results in [15]-[18]. Further we have multiphed the LO cross section by 
= 1 -I- 22 asjA'K in Eq. (2.10). Finally we have the approximation for 
arbitrary masses mu and mt 



^ ^NLO,approx 

dpt 



{S,Pt,mH) 



d a 



LO.exact 



dpt 



-{S,pl,m]j) 



J _.S-|-V,approx 

+^^^^- (5,p?,m^), (3.12) 

where the partonic cross sections are given in Eqs. (2.4), (2.6) and (2.7) 
and the soft-plus- virtual gluon approximation is given in Eqs. (2.14) and 
(2.16). The running coupling and parton densities are also represented in 
next-to-leading order for which we have chosen the MS-scheme and Uf — 5. 

For our plots we have adopted the parametrizatiou obtained from the set 
MRST [29] (namely set alfll9.dat) with A^^^^ = 239 MeV. For simplicity 
the factorization scale /x is set equal to the renormalization scale /x^- For our 
plots we take n 



1^0 



Pt 



niff unless mentioned otherwise. 
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Our first study concerns the validity of the soft-phis- virtual gluon approx- 
imation. This is done in the asymptotic top-quark limit where we know the 
complete NLO correction. For that purpose we plot 

in the range 40GeV/c < pt < 200 GeV/c and tuh = 120, 160, 200GeV/c^ 
The plots are given for the LHC (v^ = 14 GeV) in Fig. 7. The fi gure reveals 
that at rriH = 120 GeV/c^ andpt = 40 GeV/c the ratio is 1.06 and it decreases 
to about 0.9 at pt — 200GeV/c. For larger Higgs masses the ratio becomes 
closer to unity at pt > 100 GeV/c. This feature can be understood because 
at larger Higgs masses the kinematics are closer to the boundary of phase 
space. The conclusion is that in the range lOOGeV/c < pt < 200GeV/c 
we have 0.9 < R < 1.0 which indicates that the soft-plus-virtual gluon 
approximation with the prescription in [27] works rather well. This is mainly 
due to the dominance of the gg-channel and the steeply rising gluon flux 
which is even enhanced by the definition of the parton luminosity in Eqs. 
(2.20), (2.21). In the case of the TEVATRON (v^ = 2 GeV) the soft-plus- 
virtual gluon approximation works even better (see Fig. 8). In the whole 
range 40GeV/c < Pt < 200GeV/c we have 0.95 < R < 1.07. However the 
mass range is more limited i.e. uih — 120, 130, 140 GeV/c^ because at larger 
masses the cross section becomes unobservably small. This is understandable 
because at lower energies we are closer to the boundary of phase space where 
the soft-plus- virtual gluon approximation approaches the exact cross section. 
The transverse momentum distributions d a/d pt are plotted in the case of 
the LHC in Figs. 9, 10, 11 for mn = 120, 160, 200GeV/c2 respectively. The 
figures reveal the differences between the cross sections in the asymptotic rrit 
limits and the exact (approximate) cross sections. They becomes more clear 
if we plot the ratios 

da^^'^^^^/dpt ' da^^'^^^^^-^/dpt ' ^ ' ' 

which are shown in Figs. 12, 13 and 14. For the Born cross section they vary 
at Pt = 40 GeV/c from 0.93 to 1.03 for mn = 120 GeV/c^ to 200GeV/c2 
respectively. At pt = 200 GcV/c they all become about 0.8 irrespective of 
the Higgs mass. For the NLO cases these values are 0.96 to 1.28 for small 
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Pt and 0.68 to 0.75 at large pt as rriH increases from itlh = 120 GeV/c^ 
to niH = 200 GeV/c^. The exact (approximate) cross sections are always 
below those in the asymptotic limit except for ttih — 160 GeV/c^ and 
tuh — 200 GeV/c^ at small pt. There are cross over points at pt = 55 GeV/c 
and Pt = 75 GeV/c for the NLO cross sections. The NLO corrections are 
very large. This becomes clear if we look at the K-f actors defined by 



which are shown in Fig. 15. At pt — 40 GeV/c the K-factors vary from 
1.66 to 2.02 as increases from mn — 120 GcV/c^ to rriH = 200 GeV/c^ 
respectively. At larger pt values the K-factors decrease and at pt = 120 GeV/c 
they stabilise around the values 1.5, 1.55, 1.65 for ttih = 120, 160, 200 GeV/c^ 
respectively. For mn — 120 GeV/c^ the difference between the asymptotic mt 
limit and the soft-plus-virtual gluon approximation in NLO is of the same 
order as the K-factors, namely 1.5. In Fig. 16 the transverse momentum 
distributions are shown for the TEVATRON at itlh = 120 GeV/c^ and in 
Fig. 17 the ratios in Eq. (3.14) are plotted. Here the discrepancies in 
NLO are even larger. There is very little difference between small pt and 
large pt and the approximate cross section is about 0.5 to 0.8 times smaller 
than the one in the asymptotic mt limit. The Born approximations vary 
from 0.8 to 1.25 when pt ranges from pt = 40 GeV/c to pt = 200 GeV/c. 
Notice that for pt > 135 GeV/c the approximate cross section becomes even 
a little bit larger than the one in the case of the asymptotic nit hmit. The 
K-factors (see Fig. 18) are a little bit smaller than in the case of the LHC. 
At Pt = 40 GeV/c they vary between 1.5 and 1.8 and at larger pt (say 
Pt > 120 GeV/c) they are in the range 1.3 < K < 1.4. Here the discrepancy 
between the asymptotic limit and the approximate cross section in NLO is 
even larger than the corresponding K-factor. Note that the peaks in Figs. 16 
and 17 reflect the thresholds in the partonic channels described in [3]-[7]. The 
dependence of the exact Born and the soft-plus- virtual gluon approximation 
cross sections on the factorization scale is studied for m/^ = 120 GeV/c^ at 
Pt = 100, 150, 200 GeV/c. The dependence can be expressed by the following 
quantity 



(3.15) 




(3.16) 
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with fiQ = ^Jpf + m^. This quantity is plotted in the range 0.1 /iq < /i < 
10 /xo for LO and NLO in Fig. 19 for the LHC and in Fig. 20 for the 
TEVATRON both at mn = 120 GeV/c^. The LO cross sections have the 
larger values for small /i/iiq- What is very striking is the improvement in 
scale variation while going from LO to NLO. In LO there is steep behaviour 
at small fi/fio which is flattened out in NLO. At large /x//io the difference 
between LO and NLO is not so big, but still the NLO curves are flatter than 
the LO ones. Basically the same curves are also found at larger Higgs masses 
so that there is no need to show them. Finally there is a small dependence 
of N{pt, n/ Ho) on the transverse momenta in both LO and in NLO. 

Concluding our flndings we observe that the soft-plus-virtual gluon ap- 
proximation gives a good description of the exact NLO cross section (within 
90 %), when tested with — > oo cross sections. The difference between 
the asymptotic rrit limit and the soft-plus-virtual gluon approximation is 
larger than the K-factor in the case of the TEVATRON but smaller than the 
K-factor in the case of the LHC. Also the validity of asymptotic rrit limit 
depends more on the value of the transverse momentum than on the magni- 
tude of the Higgs mass. Finally our approximation has a significantly smaller 
scale dependence for both colliders in particular at small factorization scale. 
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Figure Captions 

Fig. 7. The quality of the soft- plus- virtual gluon approximation represented 
by the ratio R in Eq. (3.13) for 40 GeV/c < Pt < 200 GeV/c at the LHC 
{VS = 14 TeV) for = 120GeV/c2 (solid line), niH = leOGeV/c^ 
(dashed line), niH = 200GeV/c2 (dotted line). 

Fig. 8. The quality of the soft-plus- virtual gluon approximation represented 
by the ratio R in Eq. (3.13) for 40GeV/c < Pt < 200GeV/c at the 
TEVATRON {VS = 2 TeV) for tuh = 120GeV/c2 (solid line), tuh = 
130GeV/c2 (dashed line), uih = 140GeV/c2 (dotted line). 

Fig. 9. Differential cross sections at the LHC {y/S = 14 TeV) with niH = 
120GeV/c^. The Born cross sections da^'^'^^^^^ / dpt (dotted line) and 
da^'^'^*^°° / dpt (dot-dashed hne). Also shown are the NLO contribu- 
tions d(7^LO'^PP™^/dpt (solid line) and da^^^'"''-'°° /dpt (dashed hne). 

Fig. 10. Differential cross sections at the LHG {\/S — 14 TeV) with mn — 
160GeV/c^. The Born cross sections da^^'^"^^ / dpt (dotted line) and 

da^^'"^^^^ / dpt (dot-dashed line). Also shown are the NLO contribu- 
tions (i(7NLO''^PP'"°^/cipt (solid line) and da^^^'"^^^"^ /dpt (dashed hne). 

Fig. 11. Differential cross sections at the LHC (v^ = 14 TeV) with niH — 
200GeV/c^. The Born cross sections da^'^'^^^/dpt (dotted line) and 
da^'~''^^^'^ /dpt (dot-dashed line). Also shown are the NLO contribu- 
tions dcT^LO'^PP^^/cipt (sohd line) and da^^^^"^'^^ /dpt (dashed hne). 

Fig. 12. The factors H^^ (dashed hne) and F^^^ (solid hne) in Eq. (3.14) 
at the LHC {\fS = 14 TeV) with tuh = 120GeV/c^ 

Fig. 13. Same as in Fig. 12 for tuh = IGOGeV/c^. 

Fig. 14. Same as in Fig. 12 for tuh = 200 GeV/c^ 

Fig. 15. The K factor (Eq. (3.15)) for the LHC (v^ = 14 TeV). tuh = 

120 GeV/c^ (solid line), tuh = 160 GeV/c^ (dashed line), tuh = 200 GeV/c^ 
(dotted line). 

Fig. 16. Same as in Fig. 9 but then for the TEVATRON {VS = 2 TeV) 
and niH = 120 GeV/c^. 
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Fig. 17. Same as in Fig. 12 but then for the TEVATRON {\/S = 2 TeV) 
and niH = 120GeV/c2. 

Fig. 18. The K factor (Eq. (3.15)) for the TEVATRON (v^ = 2 TeV). 
niH = 120GeV/c^ (sohd hne), ttih = 130GeV/c^ (dashed hne), ttih = 
140GeV/c2 (dotted hne). 

Fig. 19. The scale dependence represented by N{pt, fi/ fi^) in Eq. (3.16) 
for the LHG {VS = 14 TeV) and uih = 120GeV/c2. The results are 
plotted in the range 0.1 < n/no < 10 with /Xq = + for pt = 
lOOGeV/c (solid line), pt = 150GeV/c (dashed line) pt = 200GeV/c 
(dotted line). The upper three curves are for da^'^'^^^^^ / dpt whereas 
the lower three curves are for da^^'^'^'^'^^°^ / dpt- 

Fig. 20. The scale dependence represented by N{pt, ji/ iiq) in Eq. (3.16) for 

the TEVATRON {s/S = 2 TeV) and = 120 GeV/c^. The results 
are plotted in the range 0.1 < /x//io < 10 with ^Uq = m^j + pi for pt = 
lOOGeV/c (solid line), pt = 150GeV/c (dashed line) pt = 200 GeV/c 
(dotted hne). The upper three curves are for da^'^'^^^^ / dpt whereas 
the lower three curves are for da^^'^'^^'^^°^ / dpt- 
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